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ABSTRACT 

Observations conducted with the Fine Guidance Sensor on Hubble Space Telescope (HST) providing 
high cadence and precision time-series photometry were obtained over 10 consecutive days in December 
2008 on the host star of the transiting exoplanet HD 17156b. During this time l.OxlO 12 photons 
(corrected for detector deadtime) were collected in which a noise level of 163 parts per million per 
30 second sum resulted, thus providing excellent sensitivity to detection of the analog of the solar 
5-minute p-mode oscillations. For HD 17156 robust detection of p-modes supports determination 
of the stellar mean density of (p*) = 0.5301 ± 0.0044gcm -3 from a detailed fit to the observed 
frequencies of modes of degree I = 0, 1, and 2. This is the first star for which direct determination of 
(p*) has been possible using both asteroseismology and detailed analysis of a transiting planet light 
curve. Using the density constraint from asteroseismology, and stellar evolution modeling results in 
M* = 1.285 ± 0.026 M Q , i?„ = 1.507 ± 0.012R Q , and a stellar age of 3.2 ± 0.3 Gyr. 
Subject headings: planetary systems — stars: oscillations — stars: individual (HD 17156) — tech- 
niques: photometric 



1. INTRODUCTION 

General knowledge of stellar structure and evolution 
for normal main-sequence stars has long been a mature 
topic, with forefront research driven by increasingly pre- 
cise observational constraints, and theoretical investiga- 
tions including effects such as heavy element diffusion 
and overshoot of material at convection zone bound- 
aries. Transiting planets, the first o f which HD 209458 
was d iscovered only a d ecade ago (jCharbonneau et al~l 
120001 iHenrv et all 120001 ) have provided a wealth of de- 
tailed physical information about dozens of planets. Of 
more relevance to the current paper, high quality tran- 
sit light curves also return a direct determination of the 
mean density of the host star, (p*), independent of stellar 
evolution models (e.g. iSeaeer fc Mallen-Omelas 1 120031 : 
iSozetti et aLll2007L and iWinn 1 120081 ) . Asteroseismology 
has long promised fundamentally new and precise con- 
straints on stellar structure and evolution (e.g. see re- 
view of iBrown fc Gillilandl [l99l . and with the use of 
state-of-the-art ground-based spectroscopic radial veloc- 
ities prompted largely by the burgeoning efforts for RV 
planet discovery has delivered several successes in recent 
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years (e.g. see review of lBedding fc K jcldsen 2008). The 
advent of space-based photometric missions promises to 
provide robust results from oscillations on a much larger 
number of stars, as hi nted at by the early successes from 
CoRoT (Michel 2008) . Detection of several low-angular- 
degree p-modes also provides a direct constraint on (p*) 
as will be discussed at length for HD 17156 herein. 

In this paper we present photometry of the host for the 
unusually interesting exoplanet HD 17156b given its long 
orbita l period of 21.2 days discovered by iFischer et al.l 
(I2lffl7h in a Dopp l er sur vey, and shown to have transits 
by lBarbieri et al.l ((2007). The nearly 10 days of high pre- 
cision HST photometry suffices to obtain secure detection 
of some 8 individual p-modes, an accurate determination 
of the asteroseismic large separation, and hence accu- 
rate mea sure of the stellar mea n density. A companion 
paper by iNutzman et al.1 (j2010f ) will present analyses of 
complementary HST observations through three separate 
transits of HD 17156b to fix (p„) from transit light curve 
analysis. This provides the first instance of obtaining 
such precise measures of this fundamental stellar param- 
eter from two entirely different techniques, thus provid- 
ing not only the intrinsically interesting physical mea- 
surement and associated interpretations, but also pro- 
viding a test of the two methods. We find a gratifying 
level of consistency from independent application of the 
two techniques, thus providing enhanced confidence that 
both perform as expected. 
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Selection of HD 17156 as the target and the obser- 
vations obtained for the asteroseismology part of this 
project are discussed in §2. Extensive, and unique, proce- 
dures invoked for the analysis of HS T Fine Guidance Sen- 
sor time-series photometry are discussed in §3. The evi- 
dence for, and detailed results of power spectral and re- 
lated analyses for stellar oscillations are covered in §4. §5 
presents interpretation using stellar evolution and eigen- 
frequency analyses. Comparisons with the transit light 
curve analyses are provided in §6 along with a look to 
the future prospects for similar results on a much larger 
number of systems from the Kepler Mission. 

2. TARGET CHOICE AND OBSERVATIONS 

2.1. Target Choice and Background on FGS Use for 
Photometry 

The underlying motivation of this study, anticipated 
for years in application with the Kepler Mission, was 
to simultaneously challenge the two relatively new tech- 
niques of transit light curve analysis and asteroseismol- 
ogy by comparing results for (/?*) on the same object, 
while at the same time providing a new benchmark of 
accuracy for stellar and planetary quantities. Applica- 
tion of asteroseismology, with well-posed observations, 
promises determinations of stellar density to <1%, and 
ages to <10% - very desirable constraints on stellar 
structure and evolution. However, to reach such results 
requires a large allocation of observing time, both to 
reach the necessary precisions of better than a meter per 
second if through radial velocities, or near 1 part-per- 
million (ppm) if through broadband photometry, and to 
reach the necessary frequency resolution of order 1 /iHz. 
For photometry these constraints essentially come down 
to requiring that of order 10 12 photons be collected over a 
period spanning about 10 days, and for which near Pois- 
son limited results can be maintained. The best ground- 
based attempt with photometry to date involved using a 
longitude-distributed network of 4-m telesc opes for a one 
week period in 1992 (jGilliland et al.lfl993| ) to find oscil- 
lations in a cluster of slightly evolved stars expected to 
have favorable amplitudes of a few times solar; although 
a technical success this study failed to detect oscillations. 
The capability of HST to provide successful asteroseis- 
mology of solar-like stars has long been expected to hold, 
but the necessity of dedicating 10 days of observations to 
one bright star has not supported successful applications. 
In the fall of 2008 HST had lost the use of its primary 
instruments, ACS, NICMOS and STIS, leaving only the 
imager WFPC2 with nearly 15 year on-orbit in service. 
In addition, Side A of the HST electronics communi- 
cation package required to utilize WFPC2 failed, thus 
precipitating a delay in the much anticipated recent ser- 
vicing mission to upgrade Hubble's instruments. When 
use of Side B did not initially succeed, we submitted 
this program to allow good use of the remaining science 
capa bility of HST at th at time: the Fine Guidance Sen- 
sors ([Nelan et al.11200 8'). The FGS can provide exquisite 
interferometric position determinations, and time series 
photometry through summing the counts of the 4 photo- 
multipliers (PMTs) on each FGS in POSition mode, as 
well as high angular resolution, narrow fiel d interferom- 
etry w ith TRANSfer mode observations, see lNelan et al.l 
(2008) for in-depth discussion of general FGS capabili- 



ties. 

Our application for DD time on HST was granted as 
GO/DD-11945, "Asteroseismology of Extrasolar Planet 
Host Stars" . Our photometric requirements from the 
FGS exceeded those reached in its previous use, which 
usually involved single-orbit visits tim ed to coincide with 
ingress or egress of planet transits ([Wittenmver et all 
[20051: iBean |[2008h . Extensive experience has shown with 
other instruments on HST that during the first orbit, as 
the spacecraft thermally adjusts at a new orientation, 
photometric stability is much lower than in successive 
orbits. We were initially granted a test block of 8 con- 
tiguous orbits to coincide with a transit of HD 17156b. 
As we expected, after the first orbit much better stability 
followed and these obs ervations are discussed and used in 
INutzman et al.l (]2010f ). With this successful demonstra- 
tion of FGS capability we were granted the 150 orbits 
necessary to follow HD 17156 for 10 days during which it 
was in the Continuous Viewing Zone (CVZ) of HST for 
much of this time. 

The choice of HD 17156 as a target followed from a re- 
view of all 42 then known transiting planet host stars to 
see which would hold the best prospect of asteroseismol- 
ogy with a 10 day block of dedicated HST time with FGS. 
Adopting stellar param eters from Fr ederic Pont's table 
at |http: / /w ww. inscience . ch / transits (no longer active - 
the equivalent, albeit updated information may be found 
at http://exopl anet.eu) and scalings for exp ected oscilla- 
tion properties (|Kieldsen fc Bedding 1995) resulted in a 
prediction of amplitudes at 6.9 ppm for HD 17156. With 
10 days in the HST CVZ a signal-to-noise of ~ 7 was pre- 
dicted for the highest amplitude modes (and ultimately 
close to the value reached) for these FGS observations. 
Since this S/N for p-mode detections was well above that 
predicted for any other transiting planet host, and the 
system is unusually interesting given its 21.2 day period 
compared to the more common 3-5 days of Hot Jupiters 
studied to date, the target choice was clear. 

Although the discovery paper for HD 17156 appeared 
only three years ago, the system is now very well studied. 
In addition to its 21.2 day orbit, the planet has an eccen- 
tricity of 0.68. HD 17156 is one of just three transiting 
extrasolar planets for which Rossiter-McLaughlin radial 
velocity observations through transit have determined 
the relative stella r rotation and orbital plane alignm ents 
to good accuracy (|Cochranll200l INarita et 1dl2009T ) . For 
HD 17156 the stellar rotation and orbit of HD 17156b 
are well aligned. Asteroseismology holds the promise of 
also providing information on the stellar rotation angle 
relative to the plane of the sky when rotational split- 
ting of individual nonradial oscillation modes can be de- 
te cted. For HP 171 56 the Rossiter-McLaughlin study 
of INarita etafl (|2009t ) fixes vsini at 4.07 ±0.28 kms^ 1 ', 
which coupled wit h the estimate of 1.45 Rq for HD 17156 
(|Winn et a l. 2009) leads to an estimate of the stellar rota- 
tion period of 18.6±1.4 days. Since our observations span 
only half the expected rotation period, rotational split- 
ting of nonradial mode frequencies will remain under our 
frequency resolution. (In principle the full splitting for I 
— 2, m = ±2 modes could be over-sampled by a factor 
of two with our window, but our marginal signal-to-noise 
seems unlikely to support such higher order interpreta- 
tions . )^HD_J/n56 hjasbeen_w^ 

cally ([Fischer et alj|2007t lAmmler-von Eiff et al.l 12009ft , 
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via transit light curve studies that determine stellar den- 
sity, and augmented with stellar evolution comparisons 
the stellar mass and radius separately (e.g. see most 
recent such study bv lWinn et al.l 120091 ) . and has a par- 
allax from Hipparc os with a relative error of only 5% 
(jvan Leeuwenl l2007). Table 3 will include a summary of 
T c g, log(g), and [Fe/H] from the spectroscopic studies, 
as well as the primary results from this paper. Not only 
is HD 17156 the best star for these HST observations 
in a technical sense (CVZ, high count rate possible, and 
'large' predicted oscillation amplitudes), it is also near 
the absolute top of transiting extrasolar planet hosts in 
terms of intrinsic scientific interest. 

2.2. The Observations 

Observations of HD 17156 with the Fine Guidance Sen- 
sor 2 on HST were obtained for 147 contiguous HST or- 
bits spanning 9.67 days over 22 - 31 December 2008. 
Use of FGS2r, which had not previously been used for 
science observations, and was scheduled for replacement 
(and now has been replaced) in servicing mission #4 was 
selected out of extreme caution at subjecting the photo- 
multipliers to a summed count level of order 10 12 pho- 
tons as would follow for this V = 8.172 star, since some 
level of lost sensitivity is expected proportional to total 
source exposure. In the end these observations were mea- 
sured to have reduced the FGS2r sensitivity by a quite 
acceptable level of 0.2 - 0.3%, suggesting that the better 
calibrated FGSs could have been safely used. HD 17156 
is only 0.6 magnitudes fainter than the level at which 
the FGS2r digital count registers would saturate, render- 
ing any obs ervations of a brighte r target at full aperture 
impossible (|Gilliland et alj|2009f) . Since FGS2r had not 
previously been used for science observations, and since 
none of the FGSs had been carefully calibrated for de- 
tector deadtime a few additional orbits were allocated 
to this project for the purpose of calibrating the FGS2r 
deadtime, ba ckground, sensitivity and noise properties 
(see results in lGilliland et al J 12009 ). 

The Fine Guidance Sensors use photomultipliers (4 in 
each FGS) located behind first a beam splitter based on 
the two senses of linear polarization, and then behind two 
Koesters prisms that generate an interferometric signal. 
The primary use of the FGSs are to track the pointing of 
the telescope on guide stars and provide feedback into the 
pointing control system. Used as a photometer in POS 
mode the observing FGS (FGS2r in this case) acquires 
and tracks the target in Finelock, while the two guide 
FGSs keep the target stably positioned to start with. 

The FGS data is collected at 40 Hz in each of the 
four PMTs per FGS. Examined at 40 Hz the relative 
counts from the PMTs behind each Koesters prism show 
strong (of order one percent) variations with character- 
istic timescales of about a second. The sum of the two 
PMTs for each polarization remains much more constant 
on these short timescales. We used the F583W filter cov- 
ering 435 - 715 nm at full aperture which provides the 
maximum throughput available from the FGSs. 

After applying a deadtime correction averaging 11.0% 
the mean number of counts per 30 seconds summed over 
all 4 PMTs on HD 17156 is 4.96 x 10 7 . Over the 9.67 
days of observing a total of 20208 time series points were 
extracted yielding a total count level of 1.00 x 10 12 pho- 
tons for which an expected sensitivity limit of 1 ppm 



would follow if the observations are limited only by Pois- 
son statistics. 

The observations after summing over the 4 separate 
PMTs in successive blocks of 1200 samples at 40 Hz are 
shown in Fig. 1. The data were transformed to relative 
photometry by dividing each sum by the global average. 
There are several features in the data cadence to call 
attention to. The HST orbital period often results in 
gaps. During the central 6 days of these observations 
when the target was in the CVZ for HST these gaps re- 
sult from passing through the South Atlantic Anomaly 
(SAA) when high particle-event rates require the FGS 
photomultiplier high voltage (HV) to be shut down. To- 
ward the start and end of the 10 days gaps also follow 
from occultations from the Earth as the CVZ viewing is 
lost. True CVZ blocks can be seen in the 2nd through 
7th days of Fig. 1 with up to 9 consecutive HST orbits 
of unbroken data collection. A gap occurs over HJD = 
27.05 - 27.33 when for unknown reasons (particle hit in 
electronics perhaps) lock was lost and FGS2r ceased ob- 
servation of HD 17156. During this 6 hour period FGS1 
and FGS3 continued to track their guide stars, and FGS2 
observed sky, or at least a close approximation of such 
about 26 arcseconds from the star - this would prove use- 
ful in calibrating corrections for sky background. Time 
series points for HD 17156 have only been used from 
periods when data quality flags show FGS2r remained 
in fine lock. The first data point is at HJD 22.63362 
(with leading 2454800 suppressed) with the last at HJD 
32.30393 for coverage of 9.67 days during which data was 
collected. The duty cycle (fraction of total time during 
which photons were usefully collected) during this inter- 
val was 72.6% reflecting the high observing efficiencies 
possible with HST. Clearly the data structure will intro- 
duce aliases in power spectra associated with the orbital 
period of HST and the one per day cycling of CVZ peri- 
ods unbroken by the SAA. 

The goal of this project is to detect several independent 
modes of oscillation on HD 17156, the highest amplitude 
of which is expected to be < 10 ppm, and to obtain a S/N 
of several per mode in doing so. The simply extracted 
time series shown in Fig. 1 show variations peak-to-peak 
at a range of 1%, or 1000 times larger than the peak 
intrinsic variations expected from the star. Fortunately 
(and by design, we would not have attempted this other- 
wise) the stellar oscillations for HD 17156 are expected 
to have periods of 6.5 - 13 minutes, i.e. well separated 
from the bulk (but not all) of characteristic time scales 
reflected in the systematics apparent in Fig. 1. The next 
section will explain the rather extreme reduction steps 
taken to bridge the gap between simply obtained, and 
necessary, noise levels in these data. 

3. REDUCTION OF FGS2R ASTEROSEISMIC DATA 

3.1. Removal of HST Orbital Artifacts 

We have experimented with starting the data reduc- 
tions using the individual 40 Hz level counts from each 
of the four contributing PMTs. Having found no ad- 
vantages in considering the highest frequencies, we have 
settled on analyzing 30 second (1200 samples at 0.025 
second) sums. We have also explored analyzing the pairs 
of PMTs separately, and again no useful utility in doing 
so was uncovered. 
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Fig. 1. — Time series photometry for 30 second sums over the 
full ten days with two days per panel shown starting at the top and 
progressing to the bottom. The numbers within the panels show 
start values for HJD days 23 - 32 (-2454800). Each panel has the 
same 1.0% full range, which to help guide the eye is also printed 
at the upper left. See text for further details. 

The only exception to working with 30 second sums 
over all four PMTs is that we have applied the deadtime 
correction for the FGS electronics separately for each 
PMT at the 40 Hz level, and subtracted the relatively 
unimportant detec tor background count level of ~3 per 
sample per PMT (iGilliland et all 120091). The e quation 
adopted for deadtime is (see|Gilliland et al. 2009|and ref- 
erences therein): 



C't = Cm/(1-0-Cm(T d /T i )) 



(1) 



where Ct is the true number of counts within the sam- 
pling interval of 7/ = 0.025 seconds, Cm is the recorded 
number of counts in this interval, and the Tjj have val- 
ues of 210.6, 306.0, 260.3 and 286.0 nano-seconds for the 
four PMTs. In a general sense the correction for detec- 
tor deadtime is relatively unimportant in this paper, but 
essential for the analysis of transit depths since both the 
absolute count level and changes to that scale with this 
correction. 

Continuing discussion of reductions by reference to Fig. 
1, note that after each gap in the observations there is a 
rapid rise in the relative intensity. The sequence within 
each new observation (single orbits of HST normally, up 
to 9 consecutive orbits when in the Continuous Viewing 
Zone) is that the two guide FGSs turn on and establish 
"fine lock" on guide stars, and only then is the High 
Voltage (HV) turned on for the science FGS and fine 
lock guiding established on HD 17156. This results in 40 
Hz time series on HD 17156 starting at zero before HV is 
turned on, a rapid (~0.1 second) transition from zero to 
^97% of full counts at HV turn on and an exponential 
ramp up to the full count level over ~5 minutes. To 
facilitate correcting the data during the ~5 minutes of 
stabilization at the start of each observing sequence, of 
which we have 101 occurrences in our 10 days, the 30 



second data sums are always started the same number 
of 40 Hz steps after HV turn on as evidenced by rapidly 
rising counts. The offset between HV turn on and the 
start of the first 30 second sum has been chosen as 21.5 
seconds, a somewhat arbitrary offset that puts the first 
30 second sum at about 0.995 of the ultimate count rate. 
More on treating the ramp up of counts is provided later. 

A first HST orbit which included brief TRANS mode 
observations to search for companions to HD 17156, and 
the start of time series observations in POS mode has 
been discarded, and is not shown in Fig. 1. The TRANs 
mode observations show that HD 17156 appears to be 
a point source to the resolution of FGS2r, thus a stel- 
lar companion with a projected separation greater than 
about 0.015 arcsecond with a SV smaller than about 3 
can be excluded. Since HD 17156 ha s a Hipparcos dis- 
tance of ~78 pc (jvan Leeuwenl 12007?) . 0.015 arcseconds 
corresponds to ~1.1 AU. The published RV data do not 
suggest a stellar mass body with period less than a few 
years, so between FGS and RV, companions with SV < 
3 are excluded. 

A number of features are immediately obvious from 
inspection of Figure 1: (a) The data values trend down- 
ward by about 0.3% over the 10 days. This has been 
separately established to be instrumental, at least at a 
precision level of ~0.1% by observing the FGS standard 
Upgren 69 b efore and after the 1 days which showed a 
similar drop (jGilliland et al.ll2009f ). A minor degradation 
of PMT sensitivity from exposure to 10 12 photons was 
expected, (b) The HST orbital timescale is shown by suc- 
cessive vertical tick marks (95.9184 ± 0.003 minutes) to 
guide the eye (and later reductions). Within HST orbits 
the time series shows similar waveforms of typical full 
amplitude ^0.025%. These waveforms apparently evolve 
slowly over the 10 days. (Not shown - the two PMT pair 
sums show different orbital waveforms, with the varia- 
tions in one being about twice the other, but both seem 
to show equally consistent orbit-to-orbit waveforms.) (c) 
At the start of each contiguous block of data the time 
series show a ramp up of full amplitude ~0.4% following 
the HV turn on. 

The initial approach to data reductions, as had appar- 
ently worked quite well for the more limited (and simpler 
orbital waveform) data during the Nov. 7, 2008 transit 
was based on attempting decorrelations against the x,y 
pointing records from all three FGSs, and the ratios of 
counts between PMT pairs in FGS2r itself. These are 
not shown; it quickly became evident that large features 
appear in the potential decorrelation vectors that do not 
appear in these times series data over 10 days. 

We will next show plots similar to Fig. 1 that show sev- 
eral decomposition terms for the obvious artifacts. These 
will include in Fig. 2 a minor additive correction made 
for varying sky background. Fig. 3 will document a slow 
variation in time, picking up time scales that are long 
compared to the HST orbit. Fig. 4 will show an HST or- 
bital wave form that is allowed to evolve slowly in time. 
The HV ramp up in counts will be detailed in Fig. 5. 
The factors for the slow variation in time, the orbital 
waveform, and the HV ramp are solved for in an itera- 
tive, least squares procedure. Fig. 6 will show a final 
tweak that is derived after the above iterative solution. 
We next detail these corrections. These corrections are 
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applied to the raw time series such that subtraction of 
the Fig. 2 term and a point-by-point division by the four 
terms from Figures 3, 4, 5 and 6 will take the raw values 
of Fig. 1 into the corrected values of Fig. 7. 

3.1.1. Subtraction of sky and particle backgrounds 

As noted in §2 the data gap over HJD = 27.05 - 27.33 
proved useful for deriving a proxy for variations in the 
sky background. FGS1 and FGS3 both observed stars 
much fainter than our target HD 17156. In order to keep 
the bearings in the FGSs lubricated the guide FGSs were 
rotated between stars - FGS1 alternated between two 
at V = 11.50 and 11.53 (x21.5 fainter than HD 17156) 
and FGS3 alternated between V = 12.97, 13.48, 14.31 
stars (x83 - 286 fainter than HD 17156). Changes of 
background occur for two reasons: (1) variations in the 
background light, e.g. scattered from the bright Earth 
limb, and (2) changes in the charged particle flux induc- 
ing counts in the PMTs. Complete discussion develop- 
ing FGS1 and FGS3 count rate changes as a proxy for 
FGS2r need not be recounted in detail here - FGS2r is 
no longer on the telescope and was only used for science 
observations for the current program. Suffice it to note 
that minor differences in scaling for each of FGS1 and 
FGS3 to FGS2r background changes exist, and these are 
slightly different between contributions from light and 
charged particle induced counts. Also complicating this, 
one of the FGS1 stars was an obvious variable, but for 
which the intrinsic variations could be filtered out. For- 
tunately, the 6 hour gap while FGS2r was observing only 
sky provided good data for calibrating use of FGS1 and 
FGS3 as proxies, and the implied corrections as shown 
in Fig. 2 are very minor for the extremely bright HD 
17156. The subtraction for sky is performed before any 
of the more important corrections to be detailed next 
are developed in an iterative sense. Only the occasional 
brush with outer SAA regions leads to spikes that still 
remain under 0.09%, and these are well determined. We 
believe that uncertainties in the sky subtraction are quite 
unimportant. 

3.1.2. Slow variation of count rate 

The slowly varying count rate is shown in Fig. 3 and 
has been derived by running a median filter of full width 
equal to 2.2 HST orbits over the data of Fig. 1 (this is 
done only after sky subtraction and having accounted for 
the orbital wave form and HV ramp terms to be discussed 
next). This filtering should have virtually no impact on 
stellar oscillations for which we expect time scales of 5- 
15 minutes in HD 17156, while the filter is 211 minutes 
wide. 

3.1.3. HST orbit induced orbital waveform 

The most interesting and important artifact in the data 
to be dealt with arises from highly repeatable systemat- 
ics tied to the HST orbit. The source of these variations 
is treated as unknown for these purposes, and it is un- 
known. The orbital waveform derived from the data is 
shown in Fig. 4. To derive this an assumed HST orbital 
phase is generated for all the data points (the HST orbital 
period is allowed to change later as needed to produce 
optimal results, but a good value is easily derived from 
inspection). The waveform is solved for on each of 1000 
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Fig. 2. — Sky plus counts from charged particles shown as correc- 
tions relative to the HD 17156 count level. See text for how FGS1 
and FGS3 counts are used as proxies for this. The sharp, positive 
spikes result from grazing the South Atlantic Anomaly. 
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Fig. 3. — A slowly varying long-term correction factor derived 
by running a median filter with width of 0.15 days over the data. 
The 0.2 - 0.3% decline overall likely follows from lost FGS PMT 
sensitivity, and not intrinsic variation of HD 17156. 



phase points spanning an HST orbit with weights set as 
cxp(— (<5t/10.4) 2 ) where 5t is the time in seconds that the 
center of given 30 second sums is out of phase alignment 
for the orbital phase point being solved for. Essentially, 
one makes a stack of all the data points folded on the 
assumed HST orbital period and makes a weighted aver- 
age at each phase point as the ratio of the dot product 
of data and weights divided by the sum of the weights. 
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In practice this solution either does not (initially) use 
the points impacted by HV turn on, or (later) includes 
these points after correction for the HV ramp. Likewise, 
the solution assumes the slow variation shown in Fig. 3 
has first been removed. The apparent slow variation in 
time of the orbital waveform is included by forming not 
a simple weighted sum at each phase point, but rather 
a quadratic fit over the 10 days in time. Later in the 
iteration cycle this is increased to fitting a cubic in time. 
Finally, with the orbital waveform in hand, which con- 
sists of the zero point, linear, quadratic and cubic poly- 
nomial terms at each of 1000 phase points, interpolation 
is used to provide the correction at the center of each 30 
second data point. We do not claim that this nested, it- 
erative solution for the orbital waveform provides either 
a unique, or optimal correction. However, inspection of 
Figures 1 and 4 shows that the derived orbital correc- 
tion does an excellent job of matching features, and thus 
removing the ffi>T-orbit induced systematics. 
a 
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Fig. 4. — Orbital waveform component of data corrections plotted 
on the same times and scale of Figure 1. 



3.1.4. High Voltage ramp up 

The correction factor for the HV ramp period is shown 
in Fig. 5. As discussed earlier each 30 second sum is 
initiated (to ^0.01 second rms) at exactly the same offset 
following HV turn on. To first order we assume that after 
each HV turn on successive data points experience the 
same suppression of count rate. But we can do better. 
The offsets in successive accumulated sums after HV turn 
on show a strong correlation against the length of time 
the HV was off, equivalently the length of time the PMTs 
and FGS electronics were not experiencing a high photon 
flux. After removing the linear correlation with time gap 
a minor dependence on time over the 10 days is also 
clear. In practice the solution for each successive HV turn 
on point is formed as a multi-linear regression over the 



time gap, and time (ignoring only the case following the 
large 6 hour time gap which is not consistently off by this 
amount). This solution is followed for the first 135 points 
in sequences, although most of the effect is restricted 
to the first 5-8 minutes. Following all corrections the 
scatter of points near the start of observing sequences 
is no larger than the general scatter in the time series. 
The HV ramp is corrected with essentially no residual 
error, thus restoring some 12 hours of data that would 
otherwise need to be dropped since the variations in these 
time periods contain frequency components that would 
directly affect the ability to cleanly detect oscillations of 
5-15 minutes. 
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Fig. 5. — The correction factor for HV ramp up, based on a bi- 
linear regression over the preceding time gap, and overall time for 
each successive point in observing sequences. 



3.1.5. Residual orbital corrections 

After having applied the several correction factors dis- 
cussed above, it was clear from inspection that some 
obvious, albeit much smaller, residuals remained. A 
quadratic fit to each individual data segment has been 
effective in dealing with much of the residuals. Before ap- 
plying this correction the overall rms had been reduced 
to 172.22 ppm, after applying a quadratic fit over each 
data segment as shown in Fig. 6 the fully corrected data 
are shown in Fig. 7. This is the time series that will be 
used later to form power spectra and search for evidence 
of p-modes. The time series rms here is 163.86 ppm. 
The deadtime corrected count level would imply a limit 
of 141.93 ppm at 30 seconds. Based on these numbers 
the removal of artifacts has provided data within 15% of 
the Poisson limit. 

3.2. Effect on Power Spectrum Frequency Content 

With the corrected time series in hand as shown in 
Fig. 7, we are now ready to explore evidence for coher- 
ent oscillations. With an rms of 163.86 ppm over 20202 
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Fig. 6. — Final corrections derived after the iterative solution and 
application for slow-drift, orbital waveform and HV ramp changes. 
The solution follows from a simple quadratic polynomial fit to each 
data segment. 



ing power spectrum seems only to reflect noise. Based on 
published scaling relations and knowledge of the stellar 
parameters, we would expect oscillations peaking over 
1.5-2.0 mHz at perhaps 30% higher than solar ampli- 
tude, with a large separation of some 87/iHz. The er- 
rors on stellar parameters, in particular the mean stellar 
densi ty of 0.589gcm-^ , and -0.103, +0.066 1-cr errors 
from iWinn et "aTT(|200ff ) , project to a large separation of 
87.3 /iHz with a 1-ct range of 79.3-92.0 /xHz. 
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Fig. 8. — The power spectrum for the corrected data on HD 
17156 as shown in Fig. 7. The mean background noise level near 4 
ppm 2 /fiHz is evident above 2.5 /xHz. The spectrum below 1 mHz 
is likely still contaminated by systematics. 
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Fig. 7. — Same as Figure 1, but after having removed all of the 
data artifacts discussed above. Note that relative to Figure 1 the 
vertical plot scale has changed by a factor of 5 to a full range of 
0.2%. The rms of these 30 second sums (ignoring 6 points deviant 
at more than 4-cr) is 163.86 ppm. 

data points retained (dropping 6 with 4-cr deviations), 
the noise level is expected to be 1.15 ppm globally in 
amplitude, likely lower and near 1.0 ppm at frequencies 
beyond 1.5 mHz. A power spectrum is shown in Fig. 8 
over 0.5 to 4.0 mHz. Below about 1 mHz the reduction 
steps discussed above, coupled with large systematics do 
not provide useful information on any stellar variations. 
Above about 2.5 mHz inspection shows that the result- 



The upper panel of Fig. 9 shows a restricted range 
of 1.0 to 2.6 mHz for the power spectrum. A window 
function is shown as the middle panel in Fig. 9. This 
was formed by generating a time series on the observed 
cadence with very low amplitude white noise and a large 
amplitude sinusoid at the same (1.7211 mHz) frequency 
as the highest peak in HD 17156. The most prominent 
sidelobes are at multiples of the HST orbit as expected, 
with smaller features reflecting the daily changes forced 
by SAA avoidance. 

Some of our reduction steps discussed above have been 
rather drastic. In particular the orbital waveform correc- 
tion shown in Fig. 4 will suppress any real oscillation fre- 
quencies that happen to coincide with HST orbital har- 
monics (multiples of the 173.828 /iHz orbital frequency). 
To quantify this we have injected a test signal of 60 ppm 
into the raw data (Fig. 1), then performed all of the 
reductions detailed in Figures 2-6, following this with 
a new power spectrum. This is done at every frequency 
for which the power spectrum is evaluated. The ratio 
of power at the injection frequency to the input value 
then forms the response function shown in the bottom 
panel of Fig. 9. As expected, signals at harmonics of the 
HST orbit are strongly suppressed. However, the por- 
tion of frequency phase space in which signal power is 
suppressed by more than 20% (response function <0.8) 
is only 4.7%. Knowledge of the response function should 
be taken into account in searching for evidence of p-mode 
oscillations. The response function is fixed in frequency, 
while of course the sidelobes of the window function shift 
in concert with source mode frequencies. 

4. EVIDENCE FOR STELLAR OSCILLATIONS 

The asteroseismic analysis of the HD 17156 data has 
been performed using the pipeline developed at the 
Kepler Asteroseismic Science Operations Center as de- 
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Fig. 9. — The upper panel shows the power spectrum over a 
restricted range of 1.0-2.6 mHz. The adopted frequency resolution 
of 0.25 /iHz oversamples by about a factor of 5. The middle panel 
shows the window function for the power spectrum from injecting 
a strong sinusoid at the times of observations at the frequency of 
the strongest peak in the upper panel. The lower panel shows the 
response, or transfer function of the data reduction procedure used 
in removing artifacts from these data. 

scribed in IChristen sen-Dalsgaar d et al.l ((2008). A more 
extended di scussion of K e pler p ipeline analysis tools can 
be found in lHuber et al.l f|2009T ) . The outline of analysis 
steps follows: (1) Calculation of the power spectrum. (2) 
Application of a matched filter (modified approach dif- 
fering from, but with similarities to the usual comb fil- 
ter) method to determine the large frequency separation. 
(3) Calculate the folded power spectrum using the large 
frequency separation. (4) Identify relation of I mode fre- 
quencies relative to the large separation. (5) Derivation 
of individual frequencies and fit to the asymptotic rela- 
tion. (6) Calculation of the oscillation amplitudes corre- 
sponding to the radial, I = modes. 

Oscillation frequencies for low-degree, p-modes are well 
approximated by a regular series of peaks for which 
the oscillation frequencies are given by the approximate 
asymptotic relation: 

u„i « Av (n + 1/2 + e) - D l(l + 1) (2) 

where A^o = (2 J* l dr/c)~ 1 corresponds to the inverse 
of the sound travel time across the stellar diameter, and 
closely relates to the stellar mean density via: 

Av ks 135(M*/i^) 1/2 /iHz (3) 

where M* and are the stellar mass and radius in solar 
units, and the large separation for the Sun is approxi- 
mately 135 /iHz. 

In the above equations n is the radial order, and I the 
angular degree of trapped oscillation modes. Dq is sen- 
sitive to the sound speed near the stellar core, and e is 
a correction factor absorbing minor frequency dependent 
corrections sensitive to the stellar surface layers. A re- 
cent review of the theory of solar-like oscillations may be 



found in iChristensen-D alsgaard (|2004l) . 

Although individual peaks may be visible in the power 
spectrum first shown in Fig. 8, guided by an expe c tation 
given the scaling relations of IKieldsen fc Beddinsj (1995) 
for largest amplitudes of about 7 ppm near 1.65 mHz, 
these results are in a low-SNR regime in which reliance 
on existence of evenly spaced peaks in frequency must be 
an inherent part of the process of coaxing information on 
the oscillations from the data. 

A smoothed version of the power spectrum is illus- 
trated in Fig. 10 which clearly shows the excess of 
power for HD 17156. The philosophy behind this anal- 
ysis is developed in order to avoid the stochastic nature 
of the excitation and damping of individual oscillation 
modes. In order to measure the oscillation amplitude in 
a way that is independent of these effects, Kjeldsen et 
al. (2005, 2008) have suggested a method that involves 
heavily smoothing the power spectrum in order to pro- 
duce a single hump of excess power that is insensitive to 
the fact that the oscillation spectrum has discrete peaks. 
Following Kjeldsen et al. (2008) we smoothed the power 
spectrum by convolving (in power) with a Gaussian hav- 
ing a FWHM of 4 times the large separation. 
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Fig. 10. — Smoothed power spectrum for HD 17156 shown as the 
solid line with a fit to the background from frequencies over 0.5 to 
1.2 and above 2.6 mHz shown as a dashed line. The power excess 
of modes over 1.3 - 2.4 mHz in the expected domain for oscillations 
on HD 17156 is clear. 

Following IChristensen-Dalsgaard et ail ([2008) and 
IHuber et al.1 [2009) we first determine the value of the 
large separat ion that best fits the power sp e ctrum , as was 
outlined by IChristensen-Dalsgaard et all (|2007| ). The 
idea is to match to the power spectrum a series of peaks 
described such that they follow the asymptotic relation 
(Eq. (2)). To suppress the effect of the last term in that 
equation (giving rise to the so-called small frequency sep- 
aration; see below) and avoid being sensitive to the small 
deviations from the asymptotic description, the analy- 
sis is carried out for a smoothed power spectrum P{v). 
Specifically, this is obtained from the original power spec- 
trum through Gaussian smoothing with a FWHM of 
3 /iHz. We now sum the power at uniformly spaced fre- 
quencies Vk = A^(fc/2 + eo) corresponding to the leading- 
order asymptotic expression, by calculating 

2no-Afe 

F(Av,e ) = E ( 4 ) 

k— 2uq — Afc 

as a function of trial values of Az/ and eo, for a suitable 
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central radial order no and a suitable range Afc. For 
each value of Ai> we determine the maximum J-" max (A^) 
of J-(Av, eo) as a function of eo- This defines what we 
call the matched filter response, as a function of A;a Note 
that the procedure essentially determines the almost uni- 
form separation of Av/2 between the nearly degenerate 
peaks corresponding to even and odd degrees (this will 
be plotted transformed to Ais) . In our analysis we varied 
eo between 1.0 and 1.5 (typical values found from stel- 
lar models), and took Afc = 5, to include a total of 11 
peaks in the analysis. The result of using trial values for 
the large separation between 40 and 120 /iHz is shown 
in Figure 11, for n = 19.5, 20, 20.5, 21 and 21.5. The 
maximum matched filter response is found at 83.60 /iHz 
and clearly does not depend on no within the range used. 
Note that sampling the spectrum at a frequency separa- 
tion of A.v/2, and over a fixed range of radial orders, 
leads to a single peak in the response, unlike other types 
of comb analysis. 
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Fig. 11. — The matched filter-response function showing the large 
separation Auo for HD 17156. The curves show results of invidi- 
vidual foldings of the spectra using no = 19.5, 20, . . . 21.5. At 40 
fiHz the search range is 860 ± 200 uHz, while at 120 /iHz the range 
is 2580 ± 600 /iHz. The more precisely determined value from fits 
of individual frequencies to the asymptotic relation is shown at 
83.44 /iHz with the dashed line. 

The results shown in Fig. 11 yielding A^ ~ 83.6 «Hz 
are stable and robust. In particular, the same value of 
the large separation can also be found by considering 
only the first, middle or last half of the data sets, and 
it also continues to be seen if the strongest peak in the 
power spectrum at 1721.2 uHz is arbitrarily set to zero. 
Thus the 9.7 days of HST observations of HD 17156 with 
FGS2r have provided a secure detection of solar-like p- 
modes. Further quantification and use of the large sepa- 
ration follows below. 

The next stage in the analysis attempts to fix further 
details of the oscillation frequencies. A goal here is to de- 
termine the small separation, which would provide con- 



straints on the stellar age if possible, and to provide in- 
dividual frequencies for as many modes as possible. A 
necessary associated goal will be to determine mode iden- 
tifications, i.e. corresponding n and I with specific fre- 
quencies. 

Fig. 12 shows the result of folding the amplitude spec- 
trum at 83.5 uHz (taking into account the 83.44 nHz 
splitting from fits to individual mode frequencies below) 
over the range of 1.2-2. 4 mHz. The presence of strong 
contrast in this figure is simply another way of demon- 
strating that 83.5 uHz is the correct large separation, e.g. 
folds at other values would show less contrast similar to 
the fall-off in distribution of the matched-filter response 
function of Fig. 11 away from 83.5 /iHz. The expectation 
is that modes of I = 0, 1, and 2 only will be visible, with 
the modes at I = 0, and 2 nearly degenerate except for 
the small separation term of Eq. (2). This leads to the 
expectation that a single isolated peak in this figure will 
correspond to I — 1, while a doubled peak with frequen- 
cies differing by much less than the large separation will 
represent the I = and 2 modes. Guidance from theoret- 
ical models to be further discussed in §5 has led to the 
identification of modes provided in the caption to Fig. 
12 and in Table 1. 

1.2 1 r-^ 1 , . . , . 1 1 




0.2 ' '— ^ ' ' ' < ' ' ' ' 

0. 30. 60. 90. 
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Fig. 12. — The power spectrum of Fig. 9 after smoothing by 
3uHz (FWHM) and averaging over successive slices of 83.5 uHz 
over the domain of 1.2— 2. 4 mHz. The dashed lines indicate fea- 
tures which left to right correspond to modes of I = 2, 0, and 1 
respectively. 

One might also expect to have a higher signal from 
the I = modes than from the / = 2 modes in Fig. 12 
which is not the case. Kjeldsen, Bedding & Christensen- 
Dalsgaard (2008) give amplitude ratios of I = 2/0 for the 
Sun as 0.81, 0.75, and 0.67 at wavelengths of 402, 500 and 
862 nm respectively, which for the flux weighted centroid 
of these FGS observations at about 550 nm implies an 
expected ratio of 0.74. It might also be expected that 
the number of I = modes detected would exceed 1 = 2; 
again this is not the case. In fact, due to the stochastic 
nature of the excitation it would not be surprising for 
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modes with I = 2 sometimes to exceed the I = mean 
amplitudes. Of the three stars with high S/N oscillations 
discussed in Chaplin et al. (2010), KIC 3656476 has a 
l\v nearest HD 17156 and it shows equal numbers of 
detected I = and 2 modes, and in 2 of the 4 jointly 
highest pairs of these the I = 2 amplitudes are larger 
than I = 0. 

Table 1 provides the individual frequencies that follow 
from identifying the 10 highest peaks in the power spec- 
trum. These all have a signal-to-noise ratio in excess of 
4 in the power spectrum of Fig. 9. From the asymp- 
totic relation we identify 8 of the 10 highest frequencies. 
Their degree and order are shown in Table 1. The errors 
in Table 1 are estimated from simulations of the time 
series. We created a large number of stochastically ex- 
cited modes with a mode lifetime of 3 days (similar to 
the lifetime of modes in the Sun) and estimated the accu- 
racy with which we could detect the frequencies (taking 
peak amplitudes at same SNR as in the present data) in 
a time series with same sampling as the present data. A 
lifetime longer than 3 days would have resulted in smaller 
estimated frequency errors than indicated, while the er- 
rors would be larger if the mode lifetime is shorter than 
3 days. The identification of the 8 frequencies that fit 
the comb structure should not be seen as an unbiased 
frequency identification since it relies on selection of the 
modes that fit the comb structure for 1=0, 1 and 2. We 
also assume the existence of the asymptotic relation and 
the structure of stellar model frequencies to fix the iden- 
tification. Another possible approach would be to use a 
wider frequency range when identifying the modes or use 
modes with lower amplitude. We have tried to use more 
modes and the fit to the stellar models does not depend 
(within the error bars) on the exact number of peaks in- 
cluded in the analysis. In the end we decided to use the 
10 strongest peaks as a way to ensure that we only work 
with the most significant peaks and of those 10, 8 peaks 
agreed with the asymptotic relation corresponding to the 
comb power in Figure 12. 

We have carried out a weighted least-squares fit of 
the asymptotic relation, Eq. (2), to the 8 identified 
frequencies, estimating the errors in coefficients from a 
Monte Carlo simulation. The resulting coefficients are: 
Av = 83.44 ± 0.15 (jlHz, Dq = 0.90 ± 0.19 /xHz, and the 
surface term e = 1.15 ± 0.04. 

Figure 13 shows the power spectrum from Fig. 9 af- 
ter smoothing with a Gaussian of FWHM = 3 /xHz on 
which the individual mode frequencies of Table 1 are in- 
cluded. Over the frequency range 1.2 - 2.5 mHz 8 of the 
10 highest peaks are flagged as identified. 

The final step in direct analyses of the amplitude 
spectrum is to assess the amplitude per mode from a 
smoothed power spectrum, and distribution of ampli- 
tudes with frequency following the method describ ed in 
Kieldsen. Bedding fe Christensen-Da lsgaard (2008), and 
Huber et al.l (|2009). Converting the smoothed power 
spectrum (Fig. 10) to power density by multiplying by 
the effective length of the observing run followed by fit- 
ting and subtracting the background noise and then mul- 
tiplying by the large separation divided by the number of 
p-modes peaks per radial order scaled to the sensitivity of 
radial modes we may calculate the mean power per radial 
mode. The square root is then taken in order to convert 
to amplitude per oscillation mode (radial modes). The 
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Fig. 13. — A version of power spectrum as in Fig. 9 after a Gaus- 
sian smoothing of 3 fiRz FWHM has been applied. Frequencies of 
individual modes as listed in Table 1 are indicated with vertical 
bars at the listed frequency. The numerical label of 0, 1, or 2 
provides the I value. 

peak amplitude is about 7 ppm for 1 = 0, quite consistent 
with pre-observation es timates made using the pub lished 
stellar parameters and lKieldsen fc Bedding) (JT995). The 
distribution of amplitudes for radial modes in HD 17156 
is contrasted to those in the Sun in Figure 14. The ex- 
tended, non-Gaussian distribution of amplitudes, albeit 
not well determined here, is similar to those found for 
Procyon (jArentoft et al.ll2008t) . 
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Fig. 14. — The amplitude per radial mode distribution for HD 
17156 contrasted to that for the Sun. See text for discussion. 



5. STELLAR EVOLUTION MODELS AND 
INTERPRETATIONS 

5.1. ASTEC - ADIPLS Analyses 

The stellar evolution mode ls are based on ASTEC 
(jChristensen-D alsgaard 2008a]) a nd the associated eigen- 
freque ncy analysis code ADIPLS (jChristensen-Dalsgaard! 
I2008bl) . Briefly the models use t he OPAL equa- 
tion of state ( Rogers et al.l 119961 ) and opacities 



(jlglesias fc Roger's 



1996) 



rameters (jAngulo et al 



and the NACRE nuclear pa- 
I1999D . The temperature gra- 
dient in the c o nvect ion zone was computed using the 
IBohm-Vitensd (|1958D mixing-length treatment, with a 
mixing-length of qml = 2.00 pressure scale heights, 
roughly calibrated to the corresponding solar models. In 
some cases convective core overshoot was included, over a 
distance of a ov min(r cc , H p ), where r cc is the radius of the 
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convectively unstable region and H p is the pressure scale 
height at this point; the overshoot region was assumed 
to be fully mixed and adiabatically stratified. Diffusion 
and settling of helium were treated accord i ng to the sim- 
plified formulation of IMichaud fc Proffitt] (| 1993ft . It was 
assumed that the initial abundances X and Z by mass 
of hydrogen and helium were related by X — 0.77 — 3Z, 
corresponding to a galactic enrichment Ay of helium 
determined by AY = 2Z; the values of X and Z were 
characterized by the observed [Fe/H], assuming a present 
solar surface composition with ZJX — 0.0253. 

Evolution tracks have been computed to match the 
classic al obs erved parameters as gi v en in iFischer et al.l 
(|2007f) and lAmmler-von Eiff et all (12009ft . leading to 
T eff = 6082 ± 60 K and [Fe/H] =0.24 ±0.03. A grid 
of models was computed, varying the mass between 1.26 
and 1.33 M Q in steps of 0.01 M©. with composition, 
characterized by Z , corresponding to [Fe/H] varying be- 
tween 0.18 and 0.30 in steps of 0.02, and with a ov = 0, 
0.05, and 0.1. All models included diffusion and settling 
of helium. Figure [15] shows the theoretical HR diagram 
with selected evolutionary tracks. The associated Table 
2 lists the models along the evolutionary tracks at which 
the computed eigenfrequencies best match the observed 
frequencies (see §5.2). 




Fig. 15. — Theoretical HR diagram with selected evolutionary 
tracks, corresponding to the models defined in Table 2. The '+' in- 
dicate the models along the full set of evolutionary sequences mini- 
mizing the difference between the computed and obser ved frequen- 
cies. T he box is centered on the L and T e fj as given in I Winn et al.l 
(2009), with a size matchi ng the errors on these quantities (T e g 
error from Table 3, L from lWinn et al.|[2009D . 



5.2. Fitting the Observed Quantities 

We have carried out least-squares fits of the observed 
frequencies in Table 1 and the observed effective temper- 
ature to the computed grid of models. For each model in 
the grid we computed the mean square difference of the 

model frequencies ^7°^ f rom the observed frequencies 

,(°bs). 



nl 



(obs) 



(mod) s 



N 



1 ^ 

nl 



cr(v n i) 



(5) 



where N — 8 is the number of observed frequencies and 
c(^nz) is the estimated error in the frequencies. This 
was minimized along each evolution track, characterized 
by the parameters {M, Z,a ov }. We first determined the 
model, A^min, in the evolution sequence with the small- 
est value of xt- We then assumed that the best-fitting 
model for these parameters could be obtained from the 
) by scaling, 



frequencies ^™° d ' ) (M 



(6) 



and determined r by minimizing xt- The results pre- 
sented in the following are based on these resulting min- 
imal xt along the evolution tracks. According to Eq. (3) 
the value of r so obtained determines the radius of the 
best-fitting model as R = r~ 2 / 3 i?(A4 m j n ); the remain- 
ing model quantities, including the effective temperature 

T j were then determined by linear interpolation in 
radius to this value. Finally, the departure from the ob- 
served effective temperature T g bs ^ was included in the 
final 

2 



2 2 i 

X =x v + 



' rp(obs) _ rri(mod) ' 



cr(Tee) 



(7) 



where cr(T e ff) is the standard error on the effective tem- 
perature. 

The results of these fits are illustrated in Fig. [16] A key 
goal of this analysis is to determine the mean stellar den- 
sity (p*) of the star; thus panels a) and b) show xt an d 
X 2 against (p*) for all the models in the grid. It is evident 
that minimizing xt along the evolution tracks leads to 
a narrow range of mean density, with well-defined min- 
ima at fixed a ov . It is interesting that some preference is 
found for models without overshoot, although the differ- 
ence between the three cases is of limited significance. In- 
cluding also the constraint of the observed T e g (panel b) 
produces a well-defined minimum in % 2 , shifted towards 
slightly higher (p*). To quantify the location and width 
of this minimum we have fitted a parabola to those points 
that have % 2 < 3%^^, where Xmin = 2-40 is the minimum 
over all the sequences in the grid. This parabola has a 
minimum at (p*) = 0.5301 gem^ 3 , and indicates a stan- 
dard error in (p«) of 0.0031 gem -3 . Augmenting that by 
a factor of \[2 to account for further possible systematic 
errors we arrive at our final estimate: 



(p*) = 0.5301 ±0.0044 gem" 



(8) 



The error estimate in Eq. (8) is derived from the curva- 
ture of the x 2 surface. This estimate may also be cast in 
the for m of a confidence interv al in the manner described 
by, e.g., [Kallinger et all (]2010ft . This is done by comput- 
ing the relative probability exp~ x / 2 for each model, and 
then normalizing so that the sum of the probabilities 
over all models is unity. The desired confidence interval 
is then obtained by forming the marginal distribution of 
these probabilities with respect to the stellar mean den- 
sity, and determining the range of densities that contain 
the desired fraction of the total probability. In the cur- 
rent case, x 2 as a function of < p* > is very well approxi- 
mated by a parabola; as a result, the marginal probability 
distribution is almost indistinguishable from a Gaussian, 
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and the 68% confidence interval agrees closely with the 
standard error in < > of 0.0031 g cm -3 given above. 




0.530 0.535 0.540 

(p.) (g cm" 3 ) 




0.530 0.535 

(p.) (g ™1 



ered models with masses between 1.28 and 1.32 Mq and 
Z = 0.0299 ([Fe/H] = 0.24) and 0.0338 ([Fe/H]= 0.30), 
the latter case including the model minimizing \t- For 
each pair (M, Z) we selected the sequence leading to the 
smallest xl- These are the evolution tracks plotted in 
Fig. [T5] and with properties listed in Table 2. 

To illustrate the quality of the fit of the computed fre- 
quencies to the observations, Fig. [IT] compares the ob- 
served and comp uted frequencies in an echelle diagram 
(jGrec et al.lll983l see caption), for the model in the grid 
which minimizes xt- ^ is evident that the fit is excellent. 
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Fig. 16. — Results of fitting the observed frequencies and effective 
temperature to the grid of stellar models (see text for details). 
Plusses, stars and diamonds correspond to models with a ov = 
(no overshoot), 0.05, and 0.1, respectively. Panel (a) shows the 
minimum mean square deviation xi of the frequencies (cf. Eq. [5j 
along each evolution track, against the mean density (p*) of the 
corresponding models. Panel (b) similarly shows the combined 
X 2 (cf. Eq. [TJl; here the dashed curve is a parabolic fit to those 
points that have \ 2 < 7.2 (see text). Finally, panel (c) shows x 2 
against the age for the models that minimize x 2 ; the different ridges 
correspond to the different masses in the grid, the more massive 
models resulting in a lower estimate of the age. 

As illustrated in Fig. [TBb the fit also provides a con- 
straint on the stellar age, although substantially affected 
by the spread in the mass in the grid, producing a broad 
nearly flat minimum in x 2 ■ On the basis of the plot we 
estimate the age as 2.8 ± 0.6 Gyr. 

In Fig. [I5]we have indicated the models minimizing xt 
for all the parameter sets in the grid. These clearly fall in 
a tightly confined region in the diagram, corresponding 
to the strongly constrained mean density. In order fur- 
ther to illustrate the properties of the fit we have consid- 



FlG. 17. — The asteroseismic echelle diagram comparing the 
observed frequencies of Table 1 (filled symbols with lcr error bars) 
with values from the theoretical model minimizing x 2 (cf- Eq. \E\ 
open symbols); the model (model no 8 in Table 2) has a mass M = 
1.3OM , [Fe/H] = 0.30 and an age of 2.67 Gyr. The horizontal 
axis shows frequency distribution within successive 83.6 /iHz slices, 
starting at a frequency of 68.8 pHz, while the vertical axis shows 
the starting frequencies of the slices. The short horizontal lines 
indicate the frequency intervals where the response function drops 
below 50% due to orbital filtering (cf. Fig. [9j. 

As an alternative fit to the model grid, which is inde- 
pendent of the individual frequencies, we have selected, 
on each evolution track, the model that matches the large 
separation Av. Averaging the resulting mean densities 
over those models where T g differs by less than 2a from 
the observed value yields (p*) = 0.5290 ± 0.0030 gem" 3 , 
fully consistent with the more detailed fit. 

We finally note that a potential problem in fit- 
ting solar-like oscillations is the effect of the near- 
surface layers on the frequencies which is not prop- 
erly taken into account in the adiabatic modeling con- 
sidered here. This effect is well known in analysis 
of he lioseismic data (e.g., IChristensen-D alsgaar d et al.l 
1996), where it can be isolated, owing to the avail- 
ability of modes over a broad range of de gree. 
iKieldsen. Bedding fc Christensen-Dalsgaardl (|2008[ ) sug- 
gested a procedure to eliminate the effect in the analysis 
of asteroseismic data by assuming a similar functional 
form as in the solar case, but a potentially different am- 
plitude, to be determined as part of the fit; in addition, 
the procedure results in an estimate of the mean den- 
sity of the star, through scaling from a suitable reference 
model. We have applied this to the observed frequencies 
in Table 1; this resulted in insignificant changes to the 
obtained fits and stellar parameters. 
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5.3. Comparing Asteroseismic and Planet-transit Mean 
Densities with Surface Gravity 

Asteroseismology, and transit light curve modeling 
when very precise relative photometry is available, can 
both provide very precise (even accurate) determinations 
of the mean stellar density of stars. Since density is a 
simple function of stellar mass and radius, both astero- 
seismology and transit light curve modeling tightly con- 
strain allowed choices of mass and radius in individual 
cases. However, arriving at unique estimates of the stel- 
lar mass and radius separately, as for example are needed 
if one wishes to use these to provide optimal estimates 
of the mass and radius of hosted planets as needed to 
advance and challenge the theoretical study of extraso- 
lar planet formation and evolution, requires independent 
input. 

For both asteroseismology and transit light curve mod- 
eling the classical approach to fixing the stellar mass and 
radius independently, given a measured stellar density, 
relies on stellar evolution models that best match ob- 
served constraints, typically some combination of tem- 
perature, metallicity and surface gravity (luminosity in 
our case with the known parallax, although we use lumi- 
nosity only as a consistency check) , hence bringing in in- 
ferences from spectroscopy. Accurate estimates of stellar 
masses and radii then depend on correct models of both 
stellar evolution and stellar atmospheres. How good are 
these? Figure [18] illustrates one outstanding issue, that 
may perhaps exist primarily for stars in the high metallic- 
ity domain, as is the case for HD 17156, and for many ex- 
oplanet hosts given the strong correla tion of high metal- 
licity and existence of such p lanets (jFischer fe Valentil 
120051 iValenti fc Fischerl 120051 ). S pectroscopic studies 
have fixed log(q ) to be 4.29 ± 0.06 (iFischer et al.ll2007f) 
and 4.33 ± 0.05 (jAmmler-von Eiff et al.ll2009ri . with very 
consistent determinations of T e ff and [Fe/H] between the 
two independent studies as well; see Table 3 for sum- 
maries. In Fig. HU we show the range of stellar masses 
and radii that jointly meet either the asteroseismic con- 
straint for HD 17156 where M/R 3 is a constant, or the 
constraint imposed by published spectroscopic log(g) val- 
ues which requires M/R 2 to be a constant. In this dia- 
gram an ideal result would be that the two sets of joint 
M, R constraints cross at the location of the preferred 
values obtained from other considerations. Instead, we 
are left with the result that there is no agreement, the as- 
teroseismic result for radius remains nearly 3 a away from 
the log(<?) constraint near the pre ferred mass. No t e that 
the transit light curve analyses of iNutzman et all ()2010l ) 
also impose a constraint with the functional form of that 
from asteroseismology, in relative terms this would have 
a central value that falls near, but slightly above the 
preferred asteroseismic value, with error offsets about 
5 times larger than those shown by asteroseismology. 
Clearly the spectroscopic log(g) of 4.31 is well above the 
asteroseismic (or transit light curve) value of about 4.19 
for HD 17156. Perhaps in the high metallicity domain 
of interest for HD 17156, and many exoplanet hosts, the 
stronger lines result in biases for log(g). 



5.4. Calibration of Stellar Mass and Radius to Eclipsing 

Binaries 
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Fig. 18. — The upper curves show the allowed M*,i?* values 
that satisfy the asteroseismology large separation constraint (solid 
curve prefered value of (p„) = 0.5301 gem" 3 , dashed lines the 
0.0044 error offsets from this). The lower curves show the corre- 
sponding values that satisfy the spectroscopic log(g) value of 4.31, 
with solid and dashed curves preferred and 0.04 error offsets. The 
large plus symbol shows the preferred value of M* , _R* , and al- 
lowed error range res ulting from the transit light curves analyses 
(Nutzman et al. 2010) - note the excellent consistency with aster- 
oseismology. 

Figure [18] shows that radius estimates based on spec- 
troscopic log(<?) determinations are problematic at the 
10% level, whereas estimates based on transit timing and 
on asteroseismology are mutually consistent. This con- 
sistency provides a necessary, but not sufficient condi- 
tion for using the stellar mean density as measured from 
transit light curves or from asteroseismology to provide 
the truly desired stellar masses and radii. To verify that 
we can combine a mean density measurement with stel- 
lar evolution models to obtain accurate masses and radii, 
we must test our techniques against stars for which there 
are accurate and independent measures of mass and ra- 
dius. Eclipsing binaries (EBs) comprise the only large 
set of stars for which such information exists. The mass 
and radius estimates derived from light curve and radial 
velocity measurements of EBs are almost independent 
of theory, and in many cases accuracies o f better than 
3% ca n be obtained for both quantities. iTorres et al.l 
(2009) have recently published the fundamental data for 
all known EBs meetin g this s t andar d of accuracy, a list 
containing 190 stars. iBrown (2010) applied the mean- 
density analysis to 156 of these (all with masses M* 
satisfying Q.4 Mq < M» < 5M Pl ), matching Yonsei-Yale 
(YY) models (|Yi et al.ll2001[ ) to observational estimates 
of mean density, effective temperature, and metallicity. 
Comparing radii and masses from model-matching to 
those measured directly showed that matching the YY 
models generally yielded accurate results, with system- 
atic errors estimated to be smaller than 2% in radius 
and 6% in mass. Larger errors can occur for cool and 
rapidly-rotating stars, in which surface magnetic activ- 
ity is thought to interfere with convective energy trans- 
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port, leading to radii that are larger than the YY mod- 
els predict. HD 17156 is, however, a slow rotator that is 
somewhat hotter than the Sun. For this star, systematic 
errors resulting from activity should be small, certainly 
less than 1% based on the expected rotation period of 19 
days mentioned in §2, and the lack of any photometric 
variations in excess of ~0.1% intrinsic to the star in our 
FGS photometry. 

5.5. Stellar Parameters from the Inferred Mean Density 

To obtain an independent estimate of the values of the 
basic stellar parameters and their errors we have car- 
ried out a Markov Chain Monte Carlo (MCMC) fit of 
stellar models to the observed T e g and [Fe/H] and the 
asteroseismically inferred (see Brown 20Tol for de- 
tails). The models were obtained from the Yonsei-Yale 
(YY) compilation (|Yi et al.ll200l . As for the ASTEC 
models used in the asteroseismic analysis the YY models 
use OPAL opacities and equation of state and include 
diffusion and settling of helium, although with a some- 
what different formulation. Core overshoot is included, 
with a step-function dependence of oversh oot distance 
on stellar mass (see IDemarque et al.l 120041 for details), 
such that for the models relevant to the present fits a ov 
is probably generally equal to 0.1. 

The inferred parameters are listed in Table 3. The 
mass, radius and luminosity are very close to the values 
obtained from the asteroseismic analysis (see Table 2). 
The age inferred from the MCMC analysis, 3.2 ±0.3 Gyr, 
is formally consistent with the asteroseismic value. How- 
ever, the difference of 0.4 Gyr probably reflects system- 
atic differences between the ASTEC and YY evolution 
codes. This deserves further investigation. 

6. DISCUSSION 

6.1. Summary and Comparison with Previous Results 

We reported the detection of the asteroseismic large 
separation for HD 17156, and determined the mean 
stellar density resulting from this to be 0.5301 ± 
0.0044 gem -3 . The be st determination b efore these ob- 
servations comes from I Winn et a l. (2009)) who relied on 
ground-based transit light curve analysis to fix (/?») at 
the 1 a range of 0.486 - 0.655 gem" 3 , with a preferred 
value of 0.58 9 gcm~ 3 . Ou r densi ty is significantly differ- 
ent than the I Winn et al.l ([2009) value, but well within 
their original confidence range. 

6.2. Comparison with Joint HST FGS Transit Analysis 

Of more relevance to comparing our asteroseismic re- 
sults to those from transit light curve based studies are 
the results from independent HS T FGS observations ob - 
tained as a part of this program. INutzman et al.l ([20100 
find a mean stellar density of 0.522+q qis g cm" 3 . This 
transit-based determination is 1.8 a from the asteroseis- 
mic result, while the latter is only 0.4 a off from the tran- 
sit analysis based on its larger error allowances. Clearly 
of the two extreme outcomes possible from this first com- 
parison of the two independent techniques that both 
yield what may be referred to as direct determinations 
of (p*) we are firmly in the domain of mutual confirma- 
tion, rather than a potential domain of disagreement call- 
ing for further understanding of one or both techniques. 
As noted in §5.3, however, the same cannot be said for 



consistency with the spectroscopically determined log(g) 
where there is a nearly 3 a discrepancy in the sense that 
the spectroscopic log(<?) is too large. 

6.3. Future Prospects from Kepler 

In the near future the Kepler Mission may be expected 
to provide similar results in which both asteroseismology 
and transit light curve analysis will provide simultaneous 
constraints for the mean stellar density of planet host 
stars. The three previously known exoplanets in the Ke- 
pler field of view will be observed at short cadence (58.8 
seconds) throughout the mission, thus supporting aster- 
oseismology on these targets. "Only" 512 targets may 
be followed at short cadence at any time, the bulk of 
observations for over 150,000 stars with Kepler will use 
the long cadence of 29.4 minutes which suffices for detec- 
tion of planets via transits. The Kepler throughput is a 
factor of about 5.9 higher than that for the iJSTFGS2r 
used for the observations in this paper, following from 
the use of back-side illuminated CCDs on Kepler, and 
a broad bandpass of roughly 420 - 880 nm, despite the 
factor of 6.4 relative advantage in aperture area of HST 
compared to Kepler. Furthermore, Kepler should reach 
a duty cycle of nearly 100% during the month-long ob- 
serving blocks between short breaks for telemetering ac- 
cumulated data to the ground, compared to the uniquely 
high 72.6% duty cycle reached with HST for this obser- 
vation taking advantage of a CVZ passage. This results 
in a net advantage in terms of Poisson statistics limit 
of 2.27 magnitudes for Kepler observations compared to 
those in this paper. What our HST observations have 
provided for a V = 8.17 star in 10 days, should be pos- 
sible with Kepler for a V — 10.44 magnitude star in the 
same length of time. Also weighing in favor of Kepler 
asteroseismology is an expected window function with- 
out the sidelobes resulting from the orbit of HST and 
the daily passages through the South Atlantic Anomaly. 
And of course it should be routine for Kepler to devote 
much longer observing periods to targets than was the 
case for this unusually long HST observation. The initial 
target catalog for Kepler long cadence observations con- 
tains about 2500 targets brighter than the level which 
should return Poisson limited precisions per unit time 
as good or better than those discussed in this paper. 
To reach the same S/N on expected oscillations in the 
three previously known exoplanet hosts within the Ke- 
pler field of view, as for these HST observations of HD 
17156 should take about 2 months for TrES-2 at V = 
11.4, 8 days for HAT-P-7 at V = 10.5, and 20 days for 
HAT-P-11 at V = 9.6 taking into account the expected 
oscillation properties of each. 

Coming at a time of significant stress on the HST 
project these observations required the expert assistance 
of many individuals to develop and execute. Especially 
noteworthy were the skill and tireless efforts provided by 
Merle Reinhart at STScI in expertly crafting the Phase 
II program to use all available observing time in the or- 
bits allocated, and the efforts from Mike Wenz at God- 
dard Space Flight Center for shepherding the proposal 
through pre-flight reviews and monitoring engineering 
performance during execution. We thank the STScI Di- 
rector, Matt Mountain, for the generous DD time award 
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that made these results possible. We thank Matt Holman STScI. 
and Jeff Valenti for discussion. Financial support for this Facilities: 
work was provided through program GO/DD-11945 from 
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TABLE 1 

Individual frequencies (//Hz) identified for HD 17156. 



n 


I = 


I = 1 


I = 2 


13 






1258.2 ± 1.4 


17 


1516.0 ± 1.4 






19 




1721.2 ±0.9 


1759.4 ± 1.4 


21 




1885.4 ± 1.3 




22 


1932.3 ± 1.2 


1970.4 ± 1.5 




24 






2177.9 ± 1.3 



Note. — Frequencies in /iHz of individual modes identified in HD 17156. 



TABLE 2 

Stellar evolution models fitting the observed frequencies in Table 1. 



No 


M,/M Q 


Age 


Z 


X Q 


i?*/R© 


(P*> 


T e f£ 


WL© 




x 1 






(Gyr) 








(gem 3 ) 


(K) 








1 u 


1.28 


2.936 


0.0299 


0.6803 


1.505 


0.5292 


6058 


2.74 


2.46 


2.57 


2 


1.29 


2.756 


0.0299 


0.6803 


1.508 


0.5301 


6079 


2.79 


2.43 


2.43 


3 a 


1.30 


2.557 


0.0299 


0.6803 


1.509 


0.5325 


6123 


2.87 


2.59 


3.14 


1 


1.31 


2.389 


0.0299 


0.6803 


1.512 


0.5337 


6144 


2.93 


2.65 


3.85 


5 


1.32 


2.215 


0.0299 


0.6803 


1.514 


0.5357 


6176 


2.99 


2.82 


5.47 


6 


1.28 


3.069 


0.0338 


0.6687 


1.507 


0.5264 


5987 


2.62 


2.48 


4.80 


7 


1.29 


2.865 


0.0338 


0.6687 


1.510 


0.5279 


6021 


2.69 


2.40 


3.31 


8 


1.30 


2.670 


0.0338 


0.6687 


1.512 


0.5294 


6054 


2.76 


2.34 


2.51 


9 a 


1.31 


2.474 


0.0338 


0.6687 


1.514 


0.5316 


6099 


2.85 


2.52 


2.64 


10 


1.32 


2.307 


0.0338 


0.6687 


1.517 


0.5329 


6120 


2.90 


2.56 


3.03 



Note. — Models minimizing \v (cf. Eq. [5j along the evolution tracks illustrated in Fig. 1151 The smallest value of xi is obtained for 
model 8. Models are shown for two values of the initial heavy-element and hydrogen abundances Zo and Xo, corresponding to [Fe/H] = 0.24 
and 0.30. Models indicated by superscript 'a' were computed with overshoot with a ov = 0.05, the remaining models had no overshoot. 



TABLE 3 
System Parameters of HD 17156. 



Parameter 


Value 


1-cr limits 


Comment 


Mean density, (p») (gem -3 ) 


0.5301 


±0.0044 


A 


Age (Gyr) 


2.8 


±0.6 


A 


Effective Temp, T cff (K) 
Surface gravity, log(g) (cgs) 
Metallicity, [Fe/H]* 


6082 


±60 


B 


4.31 


±0.04 


B 


0.24 


±0.03 


B 


Mass, M* (Mq) 


1.285 


±0.026 


C 


Radius, R* (R©) 


1.507 


±0.012 


C 


Age (Gyr) 


3.2 


±0.3 


c 


Luminosity, L* (L©) 


2.79 


±0.14 


D 


Surface gravity, log(g) (cgs) 


4.191 


±0.004 


D 



Note. — A: Based on asteroseismic analysis of this paper. B: Spectroscopic results avera ging from Fischer ct al. (2007) and 
Ammlcr-von Eiff ct al. (2009). * For the purposes of deriving MCMC based errors l(Brown 2010) on M* and R* we have adopted a 
doubling of the [Fe/H] error to ± 0.06. C: Function of A and B parameters and use of stellar evolution models. D: Derived from B and C 
parameters. 



